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Abstract 



A search for events with large missing transverse momentum, jets, and at least two tau leptons 
has been performed using 2 fb _1 of proton-proton collision data at y/s = 7 TeV recorded with the 
ATLAS detector at the Large Hadron Collider. No excess above the Standard Model background 
expectation is observed and a 95% CL upper limit on the visible cross section for new phenomena 
is set, where the visible cross section is defined by the product of cross section, branching fraction, 
detector acceptance and event selection efficiency. A 95% CL lower limit of 32 TeV is set on the 
gauge-mediated supersymmetry breaking scale A independent of tan/3. These limits provide the 
most stringent tests to date in a large part of the considered parameter space. 
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Abstract 

A search for events with large missing transverse momentum, jets, and at least two tau leptons has been performed using 
2fb~ 1 of proton-proton collision data at y/s = 7TeV recorded with the ATLAS detector at the Large Hadron Collider. 
No excess above the Standard Model background expectation is observed and a 95% CL upper limit on the visible cross 
section for new phenomena is set, where the visible cross section is defined by the product of cross section, branching 
fraction, detector acceptance and event selection efficiency. A 95 % CL lower limit of 32 TeV is set on the gauge-mediated 
supersymmetry breaking (GMSB) scale A independent of tan j3. These limits provide the most stringent tests to date in 
a large part of the considered parameter space. 



1. Introduction 

Supersymmetry (SUSY) [lH5J introduces a symmetry 
between fermions and bosons, resulting in a SUSY part- 
ner (sparticle) for each Standard Model (SM) particle with 
identical mass and quantum numbers except a difference 
by half a unit of spin. As none of these sparticles have been 
observed, SUSY must be a broken symmetry if realised in 
nature. Assuming i?-parity conservation [6j, [7j, sparticles 
are produced in pairs. These would then decay through 
cascades involving other sparticles until the lightest SUSY 
particle (LSP) is produced, which is stable. 

Minimal gauge-mediated supersymmetry breaking 
(GMSB) [8l4l3| models can be described by six parame- 
ters: the SUSY breaking mass scale felt by the low-energy 
sector (A), the messenger mass (M mcss ), the number 
of SU(5) messengers (N5), the ratio of the vacuum 
expectation values of the two Higgs doublets (tan/3), 
the Higgs sector mixing parameter (fi) and the scale 
factor for the gravitino mass (C grav )- In this analysis A 
and tan/3 are treated as free parameters and the other 
parameters are fixed to M nlcss = 250 TeV, N$ = 3, /1 > 
and Cg ra v = 1, similar to other GMSB benchmark points 



in the literature, e.g. G2a [14| and SPS7 (15|. The C grav 
parameter determines the lifetime of the next-to-lightcst 
SUSY particle (NLSP). For C grav = 1 the NLSP decays 
promptly (ctnlsp < 0.1mm). With these parameters, the 
production of squark and/or gluino pairs is expected to 
dominate at the present Large Hadron Collider (LHC) en- 
ergy. These sparticles decay directly or through cascades 
into the NLSP, which subsequently decays to the LSP. In 
GMSB models, the LSP is the very light gravitino (G). 
Due to the gravitino's very small mass of O(keV), the 
NLSP is the only sparticle decaying into the LSP. This 
leads to multiple jets and missing transverse momentum 
(E™ ss ) in the final states. The experimental signature 



is then largely determined by the nature of the NLSP, 
which can be either the lightest stau (fi), a right handed 
slepton (£r), the lightest neutralino (%?), or a sneutrino 
(P), leading to final states containing taus, light leptons 
(£ = e, fi), photons, b-jets, or neutrinos. For N 5 — 3 the 
fi and £r NLSPs become more dominant compared to 
lower values of N 5 . At large values of tan/3, the fi is the 
NLSP for most of the parameter space, which leads to 
final states containing between two and four tau leptons. 
In the so-called CoNLSP [16) region, the mass difference 
between the fi and the £r is smaller than the tau lepton 
mass such that both sparticles decay directly into the 
LSP and are therefore NLSP. 



pmiss 

_C/rp 



This Letter reports on the search for events with large 
jets, and at least two hadronically decaying tau lep- 
tons. The analysis has been performed using 2fb~ 1 of 
proton-proton (pp) collision data at y/s = 7 TeV recorded 
with the ATLAS detector at the LHC between March 
and August 2011. Although the analysis is sensitive to 
a wide variety of models for physics beyond the Standard 
Model, the results shown here are interpreted in the con- 
text of a minimal GMSB model. The three LEP Col- 
laborations ALEPH 0, DELPHI [HI and OPAL 
studied f\ pair production, with the subsequent decay 
fi — > tG in the minimal GMSB model. The best lim- 
its are set by the OPAL Collaboration and f\ NLSPs with 
masses below 87.4 GeV are excluded. A limit on the SUSY 
breaking mass scale A of 26 TeV was set for N$ = 3, 
M mess = 250 TeV, independent of tan ,5 and the NLSP 
lifetime. The CMS Collaboration searched for new physics 
in same-sign ditau events [20] and multi-lepton events in- 
cluding ditaus 2l| using 35pb _1 of data, but the minimal 
GMSB model was not considered. A search for supersym- 
metry in final states containing at least one hadronically 
decaying tau lepton, missing transverse momentum and 
jets with the ATLAS detector is presented in another Let- 
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2. ATLAS detector 



The ATLAS detector 23] is a multi-purpose appara- 
tus with a forward-backward symmetric cylindrical geom- 
etry and nearly Ait solid angle coverage. The inner track- 
ing detector (ID) consists of a silicon pixel detector, a 
silicon strip detector and a transition radiation tracker. 
The ID is surrounded by a thin superconducting solenoid 
providing a 2 T magnetic field and by fine-granularity 
lead/liquid-argon (LAr) electromagnetic calorimeters. An 
iron/scintillating-tile calorimeter provides hadronic cover- 
age in the central rapiditj{J range. The endcap and forward 
regions are instrumented with liquid- argon calorimeters for 
both electromagnetic and hadronic measurements. An ex- 
tensive muon spectrometer system that incorporates large 
superconducting toroidal magnets surrounds the calorime- 
ters. 



3. Simulated samples 

Monte Carlo (MC) simulations are used to extrapo- 
late backgrounds from control regions (CRs) to the sig- 
nal region (SR) and to evaluate the selection efficiencies 
for the SUSY models considered. Samples of W and 
Z/j* production with accompanying jets are simulated 
with ALPGEN Q, using CTEQ6L1 [25| parton density func- 
tions (PDFs). Top quark pair production, single top pro- 
duction and diboson pair production are simulated with 
MC@NL0 p6l-|2l| and the next-to-leading order (NLO) PDF 
set CTEQ6.6 29J. Fragmentation and hadronisation are 
performed with HERWIG [3(J, using JIMMY [3l[ for the un- 
derlying event simulation and the ATLAS MC10 parame- 
ter tune 0. TAUOLA ^^ and PHOTOS [35| are used to 
model the decays of r leptons and the radiation of photons, 
respectively. The production of multi-jet events is simu- 
lated with PYTHIA 6.4.25 J36J using the AMBT1 tune [37 1 
and MRST2007 LO* [H| PDFs. For the minimal GMSB 
model considered in this analysis, the SUSY mass spec- 
tra are calculated using ISA JET 7.80 jjji]. The MC sig- 
nal samples are produced using HERWIG++ 2.4.2 [40] with 
MRST2007 LO* PDFs. NLO cross sections are calculated 
using PR0SPIN0 2.1 J4ll - l46J |. All samples are processed 
through the GEANT4-based simulation \4J§ of the ATLAS 
detector [48]. The variation of the number of pp inter- 



actions per bunch crossing (pile-up) as a function of the 
instantaneous luminosity is taken into account by model- 
ing the simulated number of overlaid minimum bias events 



1 ATLAS uses a right-handed coordinate system with its origin 
at the nominal interaction point (IP) in the centre of the detector 
and the 2-axis along the beam pipe. The x-axis points from the 
IP to the centre of the LHC ring and the y-axis points upward. 
Cylindrical coordinates (R, <j>) are used in the transverse plane, <f> 
being the azimuthal angle around the beam pipe. The pseudorapidity 
is defined in terms of the polar angle 9 as r\ = — lntan(0/2). 



according to the observed distribution of the number of 
pile-up interactions in data, with an average of ~ 6 inter- 
actions. 

4. Object reconstruction 

Jets are reconstructed using the anti-fct jet clustering 
algorithm [49| with radius parameter R = 0.4. Their 
energies are calibrated to correct for calorimeter non- 
compensation, upstream material and other effects [501 ] . 
Jets are required to have transverse momentum (pr) above 
20GeVand \q\ < 2.5. 

Muons are identified as tracks in the ID matched to track 
segments in the stand-alone muon spectrometer, while 
electrons are identified as isolated tracks with a corre- 
sponding energy deposit in the electromagnetic calorime- 
ter. The selection criteria applied to muons and "medium" 
quality electrons are described in more detail in Rcfs. 



Refs. 51] and [52|, respectively. 



The measurement of the missing transverse momentum 
two-dimensional vector p™ lss (and its magnitude _E™ 1SS ) is 
based on the transverse momenta of identified jets, elec- 
trons, muons and all calorimeter clusters with \rj\ < 4.5 



not associated to such objects [53j. For the purpose of the 
measurement of E™ lss , taus are not distinguished from jets. 

In this search, only hadronically decaying taus are con- 
sidered. The tau reconstruction is seeded from anti-fc t 
jets with px > lOGeV. An r\- and px-dependent en- 
ergy calibration to the hadronic tau energy scale is ap- 
plied. Hadronic tau identification is based on observables 
sensitive to the transverse and longitudinal shape of the 
calorimeter shower and on tracking information, combined 
in a boosted decision tree (BDT) discriminator [54] . Tran- 
sition radiation and calorimeter information is used to veto 
electrons misidentified as taus. A tau candidate must have 
Pt > 20GeV, 1 77 1 < 2.5, and one or three associated tracks 
of pt > 1 GeV with a charge sum of ±1. The efficiency 
of the BDT tau identification (the "loose" working point 
in Ref. J54|), determined using Z — > tt events, is about 
60%, independent of pr, with a jet background rejection 
factor of 20-50. 

During a part of the data-taking period, an electronics 
failure in the LAr barrel EM calorimeter created a dead 
region in the second and third layers, corresponding to ap- 
proximately 1.4 x 0.2 radians in A77 x Acf>. Electron and tau 
candidates falling in this region are discarded. A correc- 
tion to the jet energy is made using the energy depositions 
in the cells neighbouring the dead region; events having at 
least one jet for which the energy after correction is above 
30 GeV are discarded, resulting in a loss of ~ 6% of the 
data sample. 

5. Data analysis 

The analysed data sample, after applying beam, detec- 
tor and data-quality requirements, corresponds to an inte- 
grated luminosity of (2.05 ±0.08) fb" 1 dEHiJ. Candidate 
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Figure 1: The px spectrum of the leading tau candidates in data 
(points, statistical uncertainty only) and the estimated SM back- 
ground after the pre-selection of candidate events, soft multi-jet re- 
jection and the requirement of two or more taus and no light lep- 
tons. The band centered around the total SM background indicates 
the statistical uncertainty. Also shown is the expected signal from a 
typical GMSB (A = 40 TeV, tan /S = 30) sample. 



events are pre-selected by a trigger requiring a leading jet, 
i.e. the jet having the highest transverse momentum of all 
jets in the event, with px > 75 GeV, measured at the raw 
electromagnetic scale, and E™ ms > 45 GeV [57] . In the 
offline analysis, these events are required to have a recon- 
structed primary vertex with at least five tracks, a leading 
jet with p T > 130 GeV and Ef iss > 130 GeV. These re- 
quirements ensure a uniform trigger efficiency that exceeds 
98%. 

Pre-selected events are then required to have at least two 
identified tau candidates and must not contain any elec- 
tron or muon candidates with transverse momenta above 
20 GeV or 10 GeV, respectively. To suppress soft multi-jet 
events, a second jet with j>t > 30 GeV is required. The pt 
spectrum of the leading tau candidate after pre-selection of 
candidate events, soft multi-jet rejection and the require- 
ment of two or more taus and no light lcptons is shown in 

Fig.ru 

This selection rejects almost all soft multi-jet back- 
ground events. Remaining multi-jet events, where highly 
energetic jets are mis-measured, are rejected by requiring 
the azimuthal angle between the missing transverse mo- 
mentum and either of the two leading jets A0(pij? 1SB ,jet 1 2) 
to be larger than 0.4 radians. 

The SR is defined by requiring m g > 700 GeV and 
m.^ 1 + mT£ > 80 GeV, where m c g is the effective masqj 
and to^ 1 + TO^ 2 is the sum of the transverse masse^j of the 
two leading tau candidates. The m s distribution after 



2 The effective mass m e jj is calculated as the sum of Ej? lsa and 
the magnitude of the transverse momenta of the two highest-px jets 
and all selected taus. 

3 The transverse mass mx formed by E™ 1BB and the px of the tau 

lepton (t) is defined as m T = */ '2p T E™ iss (1 - cos(A0(r,p™ iss ))). 



the A(j)(p™ lss , jet 1 2 ) requirement and the m^ + m^ dis- 
tribution after the m c g requirement are shown in Fig. [2] 
After applying all the analysis requirements, 3 events are 
selected in the data. 
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Figure 2: Distributions of variables used for the signal region def- 
inition in data (points, statistical uncertainty only) and the esti- 
mated SM background after the pre-selection of candidate events, 
soft multi-jet rejection and the requirement of two or more taus and 
no light leptons. The yellow band centered around the total SM 
background indicates the statistical uncertainty. Also shown is the 
expected signal from a typical GMSB (A = 40 TeV, tan = 30) sam- 
ple. 



6. Background estimation 

The dominant backgrounds in the SR arise from top- 
pair plus single top events (here generically indicated as 
tt),W—t tv t events and Z — > tt events. While the latter 
comprises final states with two true taus, which are well de- 
scribed in the simulation, the W and tt background consist 
of events in which one real tau is correctly reconstructed 
and the other tau candidates are mis-reconstructed from 
hadronic activity in the final state. Since mis-identified 



taus are not well described in the MC, the background 
contribution from tt and W — > tv t is determined simulta- 
neously in a CR defined by inverting the m e g cut. Owing 
to the requirement on A</> and of two or more taus, this CR 
has negligible contamination from multi-jet events. More- 
over, a totally negligible contribution is expected in this 
CR from signal events. The MC overestimates the number 
of events in the CR compared to data, due to mis- modeling 
of tau misidentification probabilities. MC studies show 
that the tau misidentification probability is, to a good ap- 
proximation, independent of m e ff, so that the measured 
ratio of the data to MC event yields in the CR can be 
used to correct the MC background prediction in the SR. 
In a similar way, the multi-jet background expectation is 
computed in a multi-jet dominated CR defined by invert- 
ing the A</> and m c g cuts. In addition, _E™ ISS /m g < 0.4 
is required to increase the purity of this CR sample. The 
extrapolated contribution of this background source to the 
SR is found to be negligible. 

7. Systematic uncertainties on the background 

The theoretical uncertainty on the MC-based corrected 
extrapolation of the W and tt backgrounds from the CR 
into the SR is estimated using alternative MC samples 
obtained by varying the renormalisation and factorisa- 
tion scales, the functional form of the factorisation scale 
and the matching threshold in the parton shower process. 
An uncertainty of 14 % is estimated from this procedure. 
Moreover, an uncertainty of 23 % is associated to the nor- 
malisation factor derived in the CR. This uncertainty is 
estimated by repeating the normalisation to data inde- 
pendently for W and tt. Systematic uncertainties on the 
jet energy scale and jet energy resolution [50| are applied 
in MC to the selected jets and propagated throughout the 
analysis, including to E!f lss . The difference in the number 
of expected background events obtained with the nominal 
MC simulation after applying these changes is taken as the 
systematic uncertainty and corresponds to 18 % each. The 
effect of the tau energy scale uncertainly on the expected 
background is estimated in a similar way and amounts to 
7%. The uncertainties from the jet and tau energy scale 
are treated as fully correlated. The tau identification ef- 
ficiency uncertainties on the background depends on the 
tau identification algorithm, the kinematics of the r sam- 
ple and the number of associated tracks. The system- 
atic uncertainties associated to the tau identification and 
misidentification are found to be 2.5% and 0.5%, respec- 
tively. For the tt and W backgrounds, these uncertain- 
ties are absorbed into the normalisation. The systematic 
uncertainty associated to pile- up simulation in MC is 1 %. 
The normalization of the Z+jets and diboson backgrounds 
is affected by the uncertainty of 3.7% on the luminosity 
This results in a 0.8 % uncertainty 



55,5 



measurement 

on the total background. The contributions from the dif- 
ferent systematic uncertainties result in a total background 
systematic uncertainty of 41 %. 



In total 5.3 ± 1.3(stat) ± 2.2(sys) background events are 
expected where the first uncertainty is statistical and in- 
cludes the statistical component of the background cor- 
rection factor uncertainty and the second is systematic. 
Roughly half of the background is composed of tt events 
and the other half is evenly split into W and Z events with 
accompanying jets. 

8. Signal efficiencies and systematic uncertainties 

GMSB signal samples were generated on a grid ranging 
from A = lOTeV to A = 80TeV and from tan/3 = 2 to 
tan/3 = 50. The number of selected events decreases signif- 
icantly with increasing A due to the reduced cross section. 
The cross section drops from 100 pb for A = 15TeV to 
5.0 fb for A = 80TeV. The selection efficiency is highest 
(« 3%) for high tan/3 and lower A values, including in the 
region of the GMSB4030 point (A = 40, tan/3 = 30) which 
is near the expected limit. It drops to 0.2% in the non-fi 
NLSP regions and for high A values. This is primarily a 
consequence of the light lepton veto and the requirement 
of two hadronically decaying taus, respectively. 

The total systematic uncertainty on the signal selec- 
tion from the systematic uncertainties discussed in Sec- 
tion [7] ranges between 7.5 % and 36 % over the GMSB 
grid. The statistical uncertainty from the limited size 
of the MC signal samples is of the order of 20%, with 
variations between 7.6 % and 59 % at the edges of the 
accessible signal range. Theory uncertainties related to 
the GMSB cross section predictions are estimated through 
variations of the factorisation and renormalisation scales 
in the NLO PR0SPIN0 calculation between half and twice 
their default values, by considering variations in a s , and 
by conside ring PDF uncertainties using the CTEQ6 . 6M PDF 
error sets [58(. These uncertainties are calculated for in- 
dividual SUSY production processes and for each model 
point, leading to overall theoretical cross section uncer- 
tainties between 6.5% and 22%. Altogether this yields 
20.8 ± 3.4(stat) ± 3.6(sys) ± 3.3(theo) signal events for the 
GMSB4030 point. 

9. Results 

Based on the observation of 3 events in the SR and a 
background expectation of 5.3±1.3(stat)±2.2(sys) events, 
an upper limit of 5.9 (7.0) events observed (expected) 
is set at 95 % Confidence Level (CL) on the number of 
events from any scenario of physics beyond the SM, using 
the profile likelihood and CL S method [59j. Uncertainties 
on the background and signal expectations are treated as 
Gaussian-distributed nuisance parameters in the likelihood 
fit. This limit translates into a 95% observed (expected) 
upper limit of 2.9 fb (3.4 fb) on the visible cross section for 
new phenomena, defined by the product of cross section, 
branching fraction, acceptance and efficiency for the se- 
lections defined in Section [5] The resulting expected and 
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Figure 3: Expected and observed 95 % CL limits on the minimal 
GMSB model parameters A and tan/3. The dark grey area indicates 
the region which is theoretically excluded due to unphysical spar- 
ticle mass values. The different NLSP regions are indicate. In the 
CoNLSP region the f\ and the tji are the NLSP. Additional model 
parameters are M meS s = 250 TeV, N5 = 3, fi > and C gra v = 1. 
The previous OPAL 19] limits are also shown. 



observed 95 % CL limits on the GMSB model parameters 
A and tan/3 are shown in Fig. [31 including the lower limits 
from OPAL [19| for comparison. These limits are calcu- 
lated including all experimental and theoretical uncertain- 
ties on the background and signal expectations. Exclud- 
ing the theoretical uncertainties on the signal cross section 
from the limit calculation has a negligible effect on the lim- 
its obtained. The best exclusion is set for A = 47 TeV and 
tan/3 = 37. The results extend previous limits and values 
of A < 32 TeV are now excluded at 95% CL, independent 
of tan/?. 

10. Conclusions 

A search for events with two or more hadronically de- 
caying tau leptons, large £"™ lss and jets is performed using 
2fb~ 1 of y/s = 7 TeV pp collision data recorded with the 
ATLAS detector at the LHC. Three events are found, con- 
sistent with the expected SM background. The results are 
used to set a model-independent 95 % CL upper limit of 
5.9 events from new phenomena, corresponding to an up- 
per limit on the visible cross section of 2.9 fb. Limits on 
the model parameters are set for a minimal GMSB model. 
The limit on the SUSY breaking scale A of 32 TeV is de- 
termined, independent of tan /3. It increases up to 47 TeV 
for tan/3 = 37. These results provide the most stringent 
tests in a large part of the parameter space considered to 
date, improving the previous best limits. 
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